CERN-TH.6834/93 
LPTENS 93/7 



Rational Theories of 2d Gravity 
from the Two-Matrix Model 



J.-M. Daul 
Laboratoire de Physique Theorique 
Departement de Physique de l'Ecole Normale Superieure 
24 rue Lhomond, 75231 Paris Cedex 05, France, 

V.A. Kazakov # and I. K. Kostov * 
Theoretical Physics Division, CERN 
CH-1211 Geneva 23, Switzerland 



The correspondence claimed by M. Douglas between the multicritical regimes of the 
two-matrix model and 2d gravity coupled with (p, q) rational matter field, is worked out 
explicitly. We found the minimal (p, q) multicritical potentials U(X) and V(Y), which 
are polynomials of degree p and g, correspondingly. The loop averages W(X) and 

are shown to satisfy the Heisenberg relations {W, X} = 1 and {W, Y} = 1 and essentially 
coincide with the canonical momenta P and Q. The operators X and Y create the two 
kinds of boundaries in the (p, q) model related by the duality (p, q) <-» (q,p). Finally, we 
present a closed expression for the two two-loop correlators, and interpret its scaling limit. 

Submitted for publication in: Nuclear Physics B 



CERN-TH.6834/93 
March 1993 

^ Permanent address: Laboratoire de Physique Theorique, Departement de Physique de l'Ecole 
Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex 05, France 

* Permanent address: Service de Physique Theorique de Saclay, CE-Saclay, F-91191 Gif-sur- 
Yvette Cedex, France 



1. Introduction 



The theory of non-critical strings, or 2d gravity with central charge of the matter 
c < 1, has made great steps forward in the last years. This theory allows two very different 
formalisms: the Liouville theory [l] and the two dimensional Regge calculus based on 
matrix models , H , Q , which are in a sense complementary to each other. Spectacular 
results in calculations of critical exponents || , @ , , ||] , @ , |10|] , macroscopic observables 

|12|, ||13|| and correlation functions ||14|| , [p~5j| , ||16|| have been accomplished. One of the 



important achievements was a non-perturbative analysis of 2d gravity in the framework of 
the double scaling limit [|I7| , which found its most elegant realization in the KdV hierarchy 



of string equations [18 



The generic rational string theory can be constructed as a chain of coupled random 



matrices []18| , ||19|| . However, one can look for a more economic way to do it. In this 
paper we address the problem of finding the simplest model of (discretized) 2d gravity 
with matter implementing the whole physics of non-critical strings; or, in other words, 
finding the simplest matrix model containing all known critical regimes corresponding to 
(p, q) conformal matter. 



The first attempt to build such a minimal description was made in ref. ||20|| , where 
various matters coupled with gravity were generated as multicritical regimes of the one- 
matrix model. Since the random surface representation (Feynman graphs) of this model 
does not contain any explicit introduction of matter (like spins at the sites of the random 
lattice) this latter appears as the multicritical excitations of "pure gravity" modes. It soon 
became clear [21], |22) that the multicritical regimes in the 1-matrix model correspond 
only to the special case (q = 2,p = 2n — 1). 



The next simplest possibility is the two-matrix model first introduced in |23[ and 
solved in the large iV limit in 



As was pointed out by M.Douglas |2q] , one can get, by 
appropriately adjusting the potentials of both matrices, the critical regimes corresponding 
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to various (p, g)-rational string theories, which in turn correspond to 2d gravity interacting 
with matter fields having central charge 

c=l-6 {p - q) \ (1.1) 
pq 

Some examples of this conjecture for particular values of p and q were considered in 



2§. 



Technically, this means that in the realization of the Heisenberg algebra [P, X] = 1/N 
in the basis of orthogonal polynomials of the two-matrix model J24] , [27], the operators 



P, X are the infinite matrices containing enough diagonals (unlike those in the one-matrix 
model) to describe the differential operators of a generalized KdV hierarchy of an arbitrary 
order. The appearance of generalized KdV flows in the two-matrix model was studied in 



ref.pql. 



In this paper, we present the explicit construction of the (p, q) critical points of the two- 
matrix model. The deviation from the critical point is due to the shift of the "cosmological 
constant" from its critical value 0, while the matter remains critical. Consider for example 
the point (3,4) describing the Ising model on a random graph ||. The critical potential 
means that we are exactly at the spin ordering phase transition point, and by shifting the 
cosmological constant (the overall coupling in the matrix model) from its critical value we 
make the volume of the system big but finite. The existence of an explicit and relatively 
simple formula for the (p, q) multicritical potentials explains why one was able to find 
explicitly the critical "temperatures" in many particular cases, such as Ising models on 
various types of random graphs. 

We start from the equations of motion, relating the canonical momenta P, Q to the 
coordinates X, Y, as they were derived (but not extensively used) in [2"5f. We wrote them 



1 Strictly speaking, in the case of a non-unitary critical point p > q + 1 the coupling of the 
matrix model is not the cosmological constant. The latter is coupled with the identity operator, 
which in this case is not the operator of minimal dimension. We will discuss this point later. 



in an operator form, valid for any N but adjusted to the double scaling limit. The method 



looks quite promising for more complicated matrix models (in [29| it was applied to the 
so-called induced gauge theory). 

We extract the critical potentials for the (p, q) critical regime from the algebraic form 
of the equations of motion. In order to achieve the (p, q) critical point it is sufficient to 
take as the potentials U(X) and V(Y) two polynomials of degree p and q, respectively: 
this allows for easier, more explicit and more general formulae than in ref. |30|. For inves- 
tigating the one-loop averages, we exploit the Heisenberg algebra for [P, X] = 1/N, which 
appears in the formalism of orthogonal polynomials. In the large iV (classical, planar) 
limit, the commutator reduces to the Poisson bracket {P, X} of the corresponding "clas- 
sical" functions. The most important observation here, from our point of view, is that 
the "KdV-momentum" P is equal to the 1-loop correlator W(x) = — ^) (up to some 
part, independent from the cosmo logical constant, giving no contribution to the Poisson 
brackets). This is true even before the continuous limit, and it gives a direct physical 
meaning to the Heisenberg relation: {W, X} = 1, allowing the direct calculation of the 
macroscopic observable W(x). 



Using all this we were able to rederive in a simple way the results of [11 1, [I2| for W(a;) 
in the general (p,q) case. The unitary (m + l,m) series can be realized in two different 
ways. The first is just the general (p, q) solution with p = m + l,q = m. The second way 
is to look for criticality under the condition that the two potentials U and V are equal. 
These two realizations turn to be physically very different. 

In the non-symmetric realization the operators X and Y have different scaling dimen- 
sions and create two different kinds of boundary. Our conjecture is that they correspond 
to the Dirichlet and Neumann boundary conditions in the statistical realization of the 
(m + l,m) model. The symmetric realization can be interpreted as a m multicritical 
regime of the Ising model on a random graph. Then the operators X and Y are related 
by the ^2 symmetry in a very simple way. 



The situation is more complicated for the non-unitary cases; what we called, in the 
unitary case, the "cosmological constant" A (the overall constant in front of the matrix 
potential) couples there with an operator of lowest negative scaling dimension A_ = 
^+f=i 1131 • The genuine cosmological constant A is coupled to the unit operator and 
measures the total volume of the random graph. These two constants are related by 

. p+i- 1 

\ = A~^~ . (1.2) 

In a recent paper [[31] the general solution of dispersionless Heisenberg equations is 
given in terms of the r-function of the Whitham hierarchy. This "solution" is not very 
explicit, but it provides a general description of the vicinity of the critical point. Therefore 



we give a presentation of our results from the point of view of ref. |31[ . 

Finally, we found the general expression for all two-loop correlators and investigated 
their scaling limit, which coincided with the known ones fpTT]] , HT2"| . In the scaling limit 



we observed a remarkable factorization property, inspiring a topological field theory in- 



terpretation. It generalizes the formula obtained in |32| in the case of the one-matrix 
model. 

This paper is organized as follows: 

In section 2 we formulate the model, and using the orthogonal polynomials formalism, 
derive the general equations of motion and their planar limit. 

In section 3 we introduce the Heisenberg algebra for the matrix eigenvalue operator 
and its momentum, obtain their planar limit in terms of the Poisson brackets and establish 
the equality between the one-loop average and the Heisenberg momentum. 

In section 4 an explicit formula for the (p, q) multicritical potentials in the two-matrix 
model will be derived. 

In section 5 we demonstrate how the equations of motion work on some examples of 
unitary as well as non-unitary theories. 

In section 6, using the scaling limit of the classical Heisenberg relation (dispersionless 
KdV) we reproduce the known results for the different loop correlators. 
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In section 7 we discuss the vicinity of the critical point, in view of our interpretation 
of the loop average as the classical canonical momentum. We review the operator content 
of the theory. 

In section 8 we interpret the two kinds of boundaries in the model as the Dirichlet 
boundaries in the two alternative standard realizations of the (p, q) critical point with q — 1 
and p — 1 spins. 

In section 9 the two-loop correlator is calculated explicitly in the two-matrix model, 
and its scaling limit is found to be consistent with the known results. 

In section 10 we summarize the results for the ^2-symmetric two-matrix model de- 
scribing the unitary series (m + 1, m). 



2. Orthogonal polynomials in the two-matrix model with two arbitrary po- 
tentials 

The partition function for the two-matrix model is defined as an integral over two 
N x N hermitean matrices X and Y 

Z = J d N2 X d N \ e ^(-u(x)-v(Y)+XY) (21) 
where U and V are arbitrary polynomial (or even non-polynomial) potentials: 

U(X) = Y^X n ,V(Y) = Y^Y n . (2.2) 
/ — : ' n ^— ' n 

n=l n=l 

The linear change of variables X — > A%X + A2, Y — > A%Y + A4 has a trivial Jacobian and 
preserves the form of the potential. Therefore the problem depends on p + q — 3 relevant 
directions in the space of parameters hi, h q , g±, g p , (3 in (|2.2| ). 

As was stated in 0, the Feynman graph expansion of this model, with respect to 
the non-quadratic terms in the potentials of ( f2.1| ) , can be interpreted as the sum over 
graphs of partition functions of Ising spins placed in the vertices of these graphs. Since 
the Ising model in two dimensions describes fermions, and the sum over Feynman graphs 
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can be interpreted as the functional integral over the two dimensional metric, it can be 
concluded that we deal with fermions interacting with 2d gravity. As is well known, the 
1/N expansion allows the classification of these partition functions according to genus, so 
that the total partition function can be viewed as the functional integral of the fermionic 
string field theory with 1/N as a string coupling. So in the limit iV — > oo we are left with 
only planar 2d manifolds. 

On the other hand, the representation of 2d manifolds in terms of planar Feynman 
graphs depends on the couplings in the potentials (|2.2| ) . One can perform the discretization 
by choosing either (p 3 - or any </> n -potential, or even any mixture of vertices of different 
orders. As in the one-matrix model ||20|| , one can adjust the corresponding couplings so 



as to get different multicritical behaviours, corresponding, in the two-matrix case, to all 
possible (p, q) minimal models coupled with 2d gravity. One can say that (unlike the case 
of SOS matter fll9ll ) the (p, q) matter is produced, in the two-matrix model, as the collective 
excitations of the Z2 (Ising) degrees of freedom and those of the discretized metric of the 
world-sheet. 

Let us start with briefly reviewing the method of diagonalization due to orthogonal 
polynomials. If we integrate with respect to the angular degrees of freedom, ( |2.1] ) becomes 
an integral over the 2N eigenvalues xx, ...,xn, y\, ..., yN'- 

Z= ]Jdx i dy i A(x)A(y)e^i=^- Uixi) - V(yi)+XiVi) . (2.3) 

J i=l 

We introduce the orthogonal polynomials \n) = ( n ( x ) an d (n\ = rj n (y) by the orthonor- 
mality relation 

(m\n) = [ dx dye^- u ^- v ^ + ^Cm(x)vn(y) = S n>m . (2.4) 



We denote the matrix elements: 

d d 

Xmn = (m\x\n), Ymn = (m\y\n) , P mn = (m\ — \n), Q mn = (m\ — \n). (2.5) 



ox ay 
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The following "equations of motion" are easily obtained by doing an integration by parts: 

jjPmn = (m\U'(x)\n) - Y mn , ^Q mn = (m\V (y)\n) - X mn . (2.6) 

At this point it is useful, with a view to the later introduction of differential operators in 
the double scaling limit, to change the notations for the indices of the matrix elements: 

X k (n) = X n - k)ni Y k (n) =Y n , n - k . (2.7) 

Now n marks the position of the matrix element on the diagonal, and k is its deviation 
from it. Then the action of the operators X, Y on the orthogonal polynomial basis is 
described by the formulae 

q-l p-1 

X\n)= X k (n)\n-k)- (n\Y = ^ Y k (n)(n-k\. (2.8) 

k=-l k=-l 

To prepare for both the continuum and planar limit, introduce the continuous variables 

t = n /f3, \ = N//3, (2.9) 

the operator h whose eigenvalues are n (i.e. h\n) = t\n),t = n/(3), and the conjugate 
coordinate u, as it was prescribed in ref. ||33|| . Then n is represented by the differential 



operator n = —(l/(3)d/du> , where (3 = ^. The constant 1/(3 plays the role of the Planck 

constant in this formalism. The operators X and P are represented by the functions 

q-i (<?-i)(p-i) 
X(u,t)= e kLV X k (n), P(u,t)= ekL ° P k(n); t = n/(3. (2.10) 

k=-l k=l 

Similarly one defines the functions Y(u, t), Q(u, t) representing the operators Y and Q 

p-l (p-l)( 9 -l) 

Y(u,t) = efe " y fc(n), Q(u,t) = e k "Qk(n); t = n/j3. (2.11) 

k=-l k=l 



It is always possible to normalize the polynomials ( |2.4| ) so that 

y_i = l,Qi = t. (2.12) 



Then the coefficient X_i is determined as a function of t = A: 



X_ 1 (n)=R(t), P 1 (n) 



R(t) 



[2.13) 



where R(t) is related to the susceptibility 



i2(A)~^logZ. 



[2.14) 



Note that for any fixed highest power of the polynomial potentials (|2.2| ), we have fixed 
highest powers with respect to z = e u in ( |2.10| ) and ( p.ll|) . 
Now eqs. (|2.6|) read 



P(u,t) = U'\X(u, t 

Q{<», t) = v'(y(u,t-^)yi-x(-u + 



ld_ 

pdi' 



(2.15) 



t -1 



The compatibility of eqs. ( 2.15 ) with the condition that P{u,t) contains only positive 
powers of e w lead, in the limit iV — > oo, to a triangular system of algebraic equations for 



the coefficients and Yjt in the expansions (|2.10| ) and (|2.11| ). 

The operator equations ( |2.15| ) are exact and valid for any N. They are very useful to 
go to the planar limit or to the double scaling limit. In the planar limit, we have to drop 
all the derivatives in the arguments of the functions X, Y in ( |2.15| ), and we arrive to the 
following system of (algebraic) equations: 



(2.16) 



P(u, t) = U (X(u, t)) - Y{-w, t) 

Q(u,t) = v'(Y(u,t))-X[-u,t). 
Again, the condition that P, Q have only positive powers of e w is sufficient to find all the 
functions P, Q, X, Y. 

In the symmetric case U = V the two equations ( |2.16| ) reduce to a single equation for 
P(u, t) = Q(u, t) and X(u, t) = Y(u, t) 



P(u, t) = U'{X{u, t)) - X(-u, t), 



[2.17) 



where, choosing the symmetric normalization of the polynomials 

t 



e~ w + ...,X = v^e w + ... (2.18) 



This equation allows the coefficients and Pk to be expressed in terms of the potential 
U. Once the potential is fixed, the first coefficient R is found as a function of t = A. For 
the (p, g)-critical potential it behaves as R ~ \ 2 /p+i- 1 4- const. 

The explicit expression of X and P can be found more easily by means of the canonical 
commutation relations, to be studied in the next section. 



3. Heisenberg algebra, its classical limit and loop correlators in the planar 
limit 



It is obvious from the definitions (2.5)of the operators X, Y, P, Q that they obey the 
Heisenberg commutation relations: 



In the classical limit iV — > 00 the commutator in eq. (|3.1| ) is replaced by the Poisson 
bracket 

\px\- dPdx dPdx -1 rs2i 

Let us now consider the loop average 

^ A) = (ra) (3 ' 3) 

as a function of u and A through x = X(u, A): W = W(X(u, A), A). We will show that 
this function has the same Poisson bracket with the coordinate X(u, A) as the momentum 
P(uj, A). By definition of the average in the basis of orthogonal polynomials: 

V(*,A) = -^/du, l—=(l) =*± . (3.4) 



d\ " ' ' 2ni J x-X(u,X) \duX\x dX 
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On the other hand, 



dW 



dX 



dW 



x dX 



dW 



d\X\u, du 



a du,X \ \ dX 



{W,X} 



(3.5) 



which implies 



{W(u,X),X(u,X)}=l. 



(3.6) 



Similarly, for 



we obtain 



W(y,X) 



Tr 1 



{W(u,X),Y(u,X)} = l. 



(3.7) 



(3.f 



Equations (|3.6| ) and (^^) mean that 



W(u,X) = P(u,X) + f(x) 



(3.9) 



W(o;,A) = g(a;,A)+^) s 



(3.10) 



where the functions f(x) and ^(x) do not depend on A. Indeed, the explicit expression for 
P as a Chebyshev polynomial f[4|] , [p4| coincides, as a function of X, with the expressions 
of the Wilson loop found in [[D] , . 

Thus the classical Heisenberg algebra (|3]6|) and (|3.8| ) is now formulated in terms of 
the basic macroscopic observables, the loop averages. 



Equation ( |3.4[ ) , which can be written as 



8P_ 
~dX 



duj 



x dX 



(3.11) 



implies that there exists a "potential" S such that 



dS 
dX 



dS_ 
dX' 



(3.12) 



The differential of this function is 



dS = PdX - ujdX. 
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(3.13) 



The function S(X, t) makes sense of the action of the Hamiltonian system, as a function 
of the coordinate X and the time t = A along the classical trajectories. 

It follows from (|3.10|) , (|3.3|) and (|3.12| ) that the function S may be obtained as the 



average 

S(x,X)=/^]og(x-X)y (3.14) 
Similarly, we could start with the y variable and the action 

S(y,\) = n±log(y-Y)\ , (3.15) 



which is, in a certain sense, dual to ( |3.14[ ). Indeed, by the equations of motion ( |2.16| ) the 



canonical pair of variables (Q,Y) is obtained from (X, P) by the involution 10 — > —uj. 

The relations of the type Q3.6| ) found in this section are, of course, valid far beyond 
the two-matrix model, since the proof ( |3.4j ) and ( |3.5|) is general enough. It is certainly 
true for the one-matrix model, for the multi-matrix chain and for any matrix model where 
the orthogonal polynomial formalism is available. One can even contemplate that these 
formulae are not specific to the orthogonal polynomial representation, and may be, some 
general approach to the string theories could be based on it. 



4. Explicit (p, g)-multicritical potentials of the two-matrix model 

The first step when investigating the (p, q) critical behaviour of the two-matrix model 
is to find a set of critical couplings gk.^k- Near such a point in the space of couplings the 
size of the planar graphs tends to infinity and the corresponding (p, q) "matter" exhibits 
long-range fluctuations. There are infinitely many ways to achieve the (p, q) critical regime 
by varying the coupling constants; we choose the most economic option, which corresponds 
to potentials U and V of minimal power p and q. Even this restriction does not fix the 
multicritical potential uniquely. We take advantage of the residual freedom to fix potentials 
with all coefficients rational. 
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At the critical point A = A* one expects the following singular behaviour of P and X 
X-X*~uj q , P - P* ~ (X - X*) p/q = u p (e U} = z) (4.1) 
Inserting this in ( |2.16[ ) , we see that the singular behaviour of Q and Y can only be 

Q-Q* ~ u q = (Y-Y*) q/p (4.2) 



The singular behaviour Q4.1|) and (|4.2| ) is sufficient to fix the form of the functions ( |2.10| ) 
and ( |2.11| ) (we normalize X so that R = 1) 

X^z) = S—^-, Y*{z) = ^—^- (4.3) 



The formal solution of (|2.16|) in terms of contour integration around z = gives the (p, q) 
multicritical potentials 

1 i X^l/z) 



p,q 2m J rj- Y*(z) 

Here the contour of integration excludes all the poles except that at z = 0. 

Note that only the negative powers of z are relevant for the definition of the potential 
from the given X*(z) and Y*(z), therefore P(z) and Q(z) do not enter the final expressions 
( fOj ) • It is easy to see that the formulae ( \i .4\ ) indeed give the polynomials of the p-th and 
q-th degree, respectively. 

The unitary series of multicritical points p = m + l,q = m can be achieved also for 
equal potentials U = V . If we take X+(z) = Y*(z) = (1 — z) m / z , then the m-critical 



potential given by (4^) will describe the point (m,m + 1) for m = 3,4, ... . For example, 
the case m = 3 will give the Ising model on a 3-coordinate random lattice (with tadpole 
and self-energy subgraphs included). 

The multicritical potentials for a few lowest p and q, as well as for the symmetric case 
for a few lowest m, are presented in Appendix A. The same calculation can be carried out 

12 



in the frequently used case of the potentials containing only even powers of matrices X 
and Y. The details are presented in section 10. In Appendix A the results for the few first 
m's in the ^-symmetric case U = V are presented. 

A remarkable feature of this symmetric realization is that for every even m the cor- 
responding multicritical potential is stable (the coefficient in front of the highest power is 
positive), which provides a proper and quite economical way to non-perturbatively define 
the double scaling limit of the model. 

5. Some examples 

Should we try to find the classical functions X{u, A), Y, P, Q when given the potentials, 
we would have to solve a triangular system of algebraic equations, with degrees as high 
as p. On the other hand, expressing the potentials in terms of the coefficients of X is 
fairly easy, as it involves only linear equations. Having recognized the critical form of 
X,Y, . . .(which are polynomials, up to some 1/z term), we find critical potentials, and we 
can then try to move away from the critical point in the coupling constants space and 
see how X, Y, . . . have to be shifted : this can be done perturbatively, for any degree. In 
this way, we may access the scaling behaviour of our macroscopic observables, and check 
whether it satisfies the Heisenberg relation. 

We shall first consider the symmetric case : the two coordinates are equal, and the 
one-loop correlator is given by the classical momentum (up to some function of X). Here, 
we reproduce the (m, m + 1) unitary case, with : 



The corresponding critical behaviour appears when we shift A away from its critical value: 
this is a particular move in the space of parameters, which involves the identity operator 
of the world-sheet in this unitary will be explained later. Here, we expect the 

susceptibility to behave as (A — Ac) 1 / 771 . 




(5.1) 



13 



The simplest situation is the Ising model on a triangulated lattice (m = 3): we may 
solve the algebraic equations ( p.!6|) explicitly and extract the scaling behaviour of X. If 
we substitute the critical form of X+(z) = (1 — z) m / z in 



P(z,t) = U'(X(z,t)) -*Q>*) 



(5.2) 



then claim that U'(<p) = a + (3 (p + 7¥? 2 , an d identify different powers of z ( we know that 
P* (z) = t* z + higher powers of z ) we may identify the critical potential: 



U'(cp) = 3-3(p-(p 2 
and the critical value t* = A* = 10. Now, we try to find X for any t : 



(5.3) 



X(z,t) 



+ a + bz + cz z 



again by identifying powers of z in ( |5.2|) and using P(z, t) = + 
following system : 

= R + c 

= 3VR + 2aVR + b 
= 3- A a -a 2 - 2 bVR 

-yjR-3b-2cVR-2ab 



(5.4) 

This leads to the 



t 



(5.5) 



R 

This way, we can express c = —R, eliminate b so as to get a quadratic equation for a with 
coefficients involving R, solve it and keep the only branch that corresponds to the critical 
value a* = —3. We obtain: 



a = 2(R-l) - V4R 2 -2R + 7 
b = -3VR-2aVR 
c= -R 



(5.6) 



where R is given through 



A = t = 32 R 3 - UR 2 + 28 R- AR(AR- 1)^AR 2 - 2R + 7 (5.7) 
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Finally, we obtain (implicitly) an exact expression for X, whose scaling behaviour is found 
by letting A — > A c and z — 1 ~ (A — A c ) ^k. We find the third Chebyshev polynomial (up 
to some arbitrary scaling on X, z ) which expresses cosh(3#) in terms of cosh(#) (see the 
next section). 

Indeed, if we consider the generic scaling limit: 



X(X, z)~(X-X c )H\ 



(z-1) 



(A-A c )^ 
P(X,z) ~ (A - A c )^r 7T 



(5.8) 

(z-1) 



(x-x c y 



with £ and n polynomials of degrees m and m + 1 respectively, and if we express the 
canonical Poisson bracket in terms of these scaling functions, we find i 

m£Ti' - (to + 1)£'tt = 1. (5.9) 

If m — 3 we find that the Chebyshev polynomials do satisfy this relation. In any case, an 
obvious computation shows that hyperbolic cosines will always satisfy this scaling relation 
(compare next section). Will it always be the only polynomial solution? 

Let us consider the m = 4 case ( critical potentials, obtained in the same way are given 
in the Appendix). Again, we try to find the coefficients of successive powers of z in the 
expression of X(z,t); the equivalent to system ( [Slf) involves equations of higher degree, 
which allows for more than one solution corresponding to the expected form for X*: one 
has to choose between different branches of solutions. Among the two possible branches, 
only one exhibits a proper critical behaviour of the string susceptibility R ~ (A — Ac) 1 / 4 . 
And the corresponding coordinate function behaves, in the critical vicinity, as the fourth 
Chebyshev polynomial. Let us note that the critical point cannot be reached by increasing 
the cosmological constant from zero towards its critical value because there is a pure gravity 
(m = 2) singularity on the way (for real susceptibilities). 



This equation was also obtained in [34] 



15 



For m = 5 the algebra becomes a bit more intricate, and among the three possible 
branches of solutions, two exhibit a fifth-order critical behaviour and have to be considered. 
One of these leads to the fifth Chebyshev polynomial for the scaling behaviour of X , while 
the other solution is not reducible to that form under rescalings of X and z; it corresponds 
to 

£(T) = T 5 + 5 T 3 + 15 T. (5.10) 

If we come back to the differential equation involving £ and it and look for all its polynomial 
solutions ( degrees 5 and 6 ) we shall find Chebyshev polynomials and another solution: 

f(T) = T 5 + 5T 3 + 15 T 

(5.11) 

tt(T) =T 6 + 6T 4 + 21T 2 + 14 
Its interpretation is not yet clear. It might correspond to another scaling behaviour of the 
matrix model around the same critical point. In any case, although m is odd here, we 
again find maxima in the plot A versus R which seem to hide the singularity of fifth order; 
again, the critical point has to be reached in a very particular way. 

Having noticed the existence of this exotic solution, we may come back to the m = 4 
case, look for all possible polynomial solutions and discover that 

£(T) =T 4 + 4 

(5.12) 

7r(T) = T 5 + 5T 

but who has to be excluded as it corresponds to the wrong scaling R ~ (A — Ac) 1 / 2 . It 
seems that all possible £, rc solutions are realized as different solutions of the algebraic 
equations associated with the matrix model. 

Let us now turn to the non-symmetric case U ^ V. As we shall see in the next section, 
it contains all (p, q) critical points (the corresponding multicritical potentials are given in 
the Appendix). 

Now, the (3,4) case reveals two branches of solution: the usual Chebyshev polynomials, 

and 

X ~ T 4 - T 

(5.13) 

P ~ 4 T 3 - 3 
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which gives the proper Poisson bracket ( exotic solution, non-existing in the symmetric 
case) and the right critical exponent for R. 

We also considered the (3,5) case, and found two possible scaling behaviours, associ- 
ated with the two possible solutions for the Poisson bracket equation in the scaling regime. 
Here Chebyshev polynomials are no longer possible solutions (as will be explained later) 
but instead: 

X = - T 5 + 5 T 3 ± 5 T 2 + 15 T ± 30 

2 (5-14) 
Y = -T 3 + 6T=f6 

these being two solutions of ( |3.2| ) ( to verify this, do not forget to change signs, as P(uj) 
scales as Y(—u) ). 

A solution for P, Q, X, Y in the unitary case (p = q + 1) will be presented in the next 
section. The complete classification of scaling behaviours is far from being accomplished. 



6. Explicit solution of the Heisenberg relation near the (p, q) critical point 



Once the scaling laws are established, it is more convenient to work with the differential 
form of the equations of motion ( |3.2|) . Let us use the parametrization 



u = 2/j, cosh(6>) 



(6.1) 



Near the (p, q) multicritical point, X ~ Q ~ yfl and Y ~ P ~ . We know by eq. 
flO|) that A ~ n p+q ~ 1 . We introduce two complete systems of polynomials in \i and u> , 
homogeneous of degree p and q 



(p) 



2 / u p cosh(/c#), k = 0,l,...,p 



n ( s q) = 2// g cosh(s0), s = 0,1,..., g 
These polynomials have Poisson brackets 



(6.2) 



p+q-l 



dX/d/j 



. sinh[(s — k)9] , , . sinh[(s + k)6] 
(ps + qk) . ' 1 + (pa - qk) 



sinh^ 



sinh^ 



(6.3) 
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Thus, expanding the functions X and P as linear combinations of ( |672| 



v 

,(p) 



p = Y,A k n[ p \ X = ^BM q) (6-4) 



fc=0 s=0 



Y = J2(-) k A k ni P \ Q = J2(-) S B S ^ (6.5) 

fe=0 s=0 

we find from (|3.2| ) a linear system of equations for the coefficients. 

This method works particularly well in the unitary case q = m,p = m + 1. Then we 
have simply 

P = ^(p) = 2^ m+1 cosh[(m + 1)6], 

X = n W = cosh(m6>) , ( 6 .6) 

A = — ^ p+q - 1 = (m+ l)fi 2m 

p + q — 1 

7. The vicinity of the critical point 

The interpretation of the momentum as the Wilson loop variable is valid also outside 
the critical point. In this section we will sketch the situation when the potentials U and V 
represent an arbitrary perturbation of a critical potential. For more details see the original 
paper by Krichever |31] . 



Let us parametrize the vicinity of the critical point by the coupling constants 
t\ = X,t2,ts,... Then the operator X will be represented in the scaling limit by a gen- 
eral polynomial of degree q. As before we introduce a "potential" S(u,ti,t2, ■■■) which in 
our model has the meaning of the partition function of a disk without marked point [eq. 
([Dl])] . The function S{X,t 1} ...) can be thought of as the action of a Hamiltonian system 
depending on the co-ordinate X and the "times" t m ,n = 1,2, ... The total differential of 
the action is 

dS = PdX - H ndtn (7.1) 



where the momentum 
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is the Wilson loop average and H n are the Hamiltonians corresponding to the times t m . 
By the definition (|7.1| ) : 

H n = -B- (7.3) 



n ~ dt 

- x 

From the commutation of the second derivatives it follows that 



dP 

dt k 



° Hh k 1.2 q 1 (7.4) 



x dX 

The relation ( |3.4j ) is a particular case of ( |7.4J ) , with k = 1. Let us mention that disper- 
sionless KP flows are determined by the Hamiltonian equations of motion 

BP BX 

Each coupling constant tk is coupled with the corresponding scaling operator Ok- Equation 
( [7.5| ) implies the zero-curvature condition 

In the construction of Krichever |3^] , the starting point was the zero-curvature condition 



( [7T5| ) which is equivalent to 

d m H n \x = d n H m \X (7.7) 

and thus implies the existence of the action such that H n = d n S\x- 

Let us consider first the unitary case p—m+1, q = m. The function H n can be thought 
of as the partition function of the disk with a puncture on the world-sheet where the O n 
operator is placed. For n = 1, 2, m— 1 this is the n th order parameter of the model. The 
operator V = 0\ is the puncture operator corresponding to the identity operator in the 
conformal theory (m, m + 1). The dimension of the time t n is that of the power // 2m+1_n , 
near the point t n = 5 n ,i/U 2m . After taking the derivative in X in order to create a marked 
point on the boundary, according to ( [7.5|) one obtains the loop-point correlators that were 
calculated in [|36| for tk = fJ? m 8k,i'- 



9 H k (X,h = A,f 2 ,...) = *i - !-(^ +fc cosh[(m + £)#]) (7.f 



dX sinh(m#) d[i 
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Integrating this formula, we find for t n = n 2m 5 n ^: 

H n = p n cosh(n#) (7.9) 

In particular, Hi = uj (this is true for arbitrary t^). If we interpret the quantity (|7.4|) in 
terms of an operator O m acting on a disk with a puncture, then it should be proportional to 
A i- 7str / 2 -(i-A m ) = x i+i/2n-(i-(m-i)/2n) and it f u ows t h a t the grav itational dimension 

of the operator Ok coupled with the "time" tk is 

A* = ^ (7-10) 
2m 

Therefore, the operators Ok, k = 1,2, ...,m are the order parameters in the (m, m + 1) 
conformal field theory coupled with gravity. The operator O m +i can be identified with 
the boundary operator d/dX marking a point at the boundary of the loop. Note that 



the functions (|7.8|) have been interpreted in [12] as the wave- functions of the on-shell 
states of the closed string. They appear as leg factors in the multi-loop string amplitudes. 
In particular, the derivative of the two-loop correlator in the two-matrix model will be 
expressed as a sum of products of these functions. 

In the general (p, q) case, the coupling A = ti no longer corresponds to the puncture 
(identity) operator, but to the operator with the minimal (negative) dimension in the 
model. This leads, as we have seen in the previous section, to more complicated expressions 
for P and X. 

There exists, however, a direction in which the solution is just as simple as that in 
the unitary case. It is along the coupling constant A = t p - q corresponding to the identity 
operator O p - q . It seems to us that this is the most natural definition of the non-unitary 
critical points. It appears when the (p, q) critical points are constructed as integrable ADE 
models on a fluctuating lattice ||36|| , ||37|| . Then the operators with negative dimensions 
do not appear by construction and the cosmological constant is coupled with the volume 
of the world-sheet. 
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Along the direction % = 5k, p - q A in the space of the coupling constants, the potential 
depends explicitly on the coupling t = A and the quasi- classical limit of (|3.1|) is no longer 
given by the Poisson bracket relation (|3.2|) . Instead, the classical functions P and X satisfy 
the new Poisson bracket condition 

dPdX dXdP fn , 

where the new variables O and A are related to the old ones by a non-canonical transfor- 
mation 

= // P-9 C osh[(p-g)0],A = fi 2q (7.12) 
Equation ( |7.11| ) follows from the explicit solution for the loop average in the (p, q) 



critical point of an ADE model on a random lattice [3 c 



P(w, A) = // p cosh(p0), X(u, A) = (j, q cosh(q6); A = fi 2q , u = /xcoshfl (7.13) 



8. Boundaries and boundary operators 

In this section we will give an interpretation of the loop operators in the two-matrix 
model in terms of statistical systems on random surfaces. 

For a string physicist Tr(x — A) -1 represents an operator creating a hole in the world- 
sheet of the string. In the two-matrix model (with two different potentials) there are two 
different loop operators 

These two operators have different dimensions; their mean values read (in the unitary case 
p = q + 1) 

w(x) = (x + Vx 2 - \) p/q + (x- Vx 2 - \y /q 

(8.2) 

W(y) = (y+ vV - \p/i)) q /ri + ( y - _ X p/i)i/p 
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The experience with the ADE and SOS models on fluctuating lattices [j3"6] , | TT| , [|3Tf allows 



usto recognize in X and Y the operators that create the two possible kinds of boundary 
in the (p, q) model. It has been noticed that the ADE models exist in two phases: dense 
and dilute. Thus each (p, q) = (q, p) model has two different statistical realizations. In one 
of them the boundary has dimension 1, and in the other a non-classical dimension p/q. 
The same is true in the case of the 0(n) model J38fl . The expansion of the two kinds of 



macroscopic loops in terms of local operators contains the two kinds of boundary operators 
in the model (see the discussion in section 2.2 of ref. ||38|| ). 

Let us illustrate this by the simplest example p = 3, q = 2. If we perform the gaussian 
integration in Y, the resulting matrix model will give the standard formulation of the 
C = theory (2,3) (pure gravity). On the other hand, it is known |3^] that this matrix 
model is equivalent to the Q = 1 Potts model on a fluctuating lattice. If we choose 
<7i = hi = 0, #2 = 1, <73 = —g , then the temperature T and the "cosmological constant" (3 
are given by 

T = l/ln/i 2 , (3 = \vL{h 2 -l)-{2/Z)\ng (8.3) 



Integrating over X we find, in the large- N limit, the following one-matrix model |3£ 

r N 

z = j U d y^ e ~ pVeff{yi) U^-y^U(^^+V^+^y 1/2 

i=1 i¥=i hi (8.4) 

Veff(y) = h 2 TrY 2 + 6gTrY + ^Tr^{^YY+ ^Tr ^_L_ T?J(^Y) 

where the constant a is determined by 

a=±-^(*-Y)- 1/2 (8.5) 

This model has the same singularity for the free energy as the standard one-matrix model 
for pure gravity, but different seal ing for the loop operators. 

Another interesting problem is to find the amplitudes involving two kinds of bound- 
aries 

W(x,y) = (Tr 1 1 ) (8.6) 
x — X x — Y 

whose calculation involves the angular degrees of freedom. 



Such amplitudes have been studied recently in |40 
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9. Two-loop correlators 

The two-point functions appear in matrix formulation as expectation values of prod- 
uct of traces; -^TrX 771 corresponds to a puncture operator whose action on the string 
world-sheet leaves a hole with a boundary of length m. Boundaries of any size may be 
considered, with the operator j?Tr , which involves all possible punctures, weighted 
with the exponential of the length of their boundaries. Such an operator has well-defined 
expectation values, provided the one-dimensional cosmological constant x is large enough. 
The continuum limit is then reached when we let x decrease towards its critical value, so 
as to allow larger and larger boundaries. 

There are three types of the loop-loop correlation functions in our model: 

Wxvix, y) = {Tr—^ Tr— ^) c , (9.1) 
x — X y — Y 

W xx (x 1 ,x 2 ) = (Tr— Ur- Tr— Ur) c (9.2) 
X\ — X X2 — X 

WVy(yi, 2/2) = (Tr— ^ Tr— U;) c (9.3) 

yi-Y y 2 -Y 

They may be computed in an explicit way by means of the orthogonal polynomial tech- 
nique. Let us start with the first correlator ( |9.1|) . In the basis of orthogonal polynomials, 
we have 

i<N-l,j>N x y 1 

Only the levels located near the Fermi surface t = A contribute to this sum, so that, in 
the planar limit, we can express the former matrix elements in terms of the functions 
X(z),Y(z): for instance, (i + 5\ is the coefficient of z~ s in — ^rr- Thus we can 

V /> ^ / ' \ 1 x—X 1 ' x-X(z) 

express the connected two-loop function as a contour integral and perform the sum over 6 

= * W X -X(z)y- Y(l) 

,5 - 1 (9.5) 
1 f f dzdz 1 1 



4tt 2 J J (l-zz) 2 x-X(z)y-Y(z) 
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where \z\ < l,\z\ < 1 for the sum to converge. 

Now, provided that x and y are big enough, for the equations X(z) — x and Y(z) = y 
to have unique solutions z(x) and z(y) for \z\ < 1, \z\ < 1, we calculate the contour integrals 
by the residua and finally obtain 

W XY (x,y) = Z ' {X ^\, 2 = -d x d y \n [1 - z(x)z(y)} (9.6) 
(1 - z{x)z{y)Y 

Note that the double pole at z = 1/z does not contribute since it gives the total derivative 
under the second integral. 

In a similar way we obtain the other two correlators 

W{x u x 2 ) = -d Xl d X2 In [ *(**)- sfo) ] (9 7) 

L X\ — X2 

W( yi ,y 2 ) = -d yi d y2 In [ %l) ~ %2) ] (9.8) 

yi - 2/2 

We see from these formulae that the knowledge of z(x) and z(y) or, in other words, of 
the one-loop averages, allows the immediate reconstruction of the corresponding two-loop 
correlators. 

Formulae ( |9.6|) - work, in principle, beyond the two-matrix model and are appli- 



cable to any matrix chain where the orthogonal polynomial formalism is valid. 

In the scaling limit z = — > l,u — > we find in the unitary case p = q + 1, using 
the explicit solution ( |6.6|) , 



and 



W XY (X, Y) = -d x d Y In [cosh(#) + cosh(r)] (9.9) 

r cosh(^i) - cosh(6> 2 ) ■, min\ 

W X x(X u X 2 ) = -d Xl d X2 In — — ^ tt-A (9.10) 

L cosh(g6'i) — cosh((/t/ 2 ) 



L cosh(gri) — cosh(gr 2 ) 



where the following parametrization is used 



X = 2/j, q cosh(g0), Y = 2^ cosh(pr) (9.12) 
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Expressions (|9.9|) - (|9.11| ) remain true even in the non- unitary case p — q > 1, if the deviation 
from the critical point is along the constant A coupled with the identity operator (see the 
discussion at the end of section 7 ). 

Let us try to relate the above results with what is known of the two-loop correlators 
for strings with discrete target space. 

The loop-loop correlator in the one-matrix model with general potential has been 
calculated in |3^] (see also the Appendix of |37[| ). In all critical regimes (q = 2, p odd) it 
is given by 

W(X l ,X 2 ) = - — ^ log(v / X 1 + ^ 2 + y/X 2 + ^) (9.13) 

which is identical to the formula ( |9.10|) , which we obtained for the XX correlator with 
X = 2cosh(#/2). The derivative of this correlator with respect to A = fx p+1 gives the 
loop-loop-point correlator and factorizes to a product of two factors 

■^-W(X 1 ,X 2 ) = — 1 - (9.14) 

dfi v /x 1 + ^ v /x 2 + ^ 



It has been noticed [|37| that such a factorization is a general property of all strings with 



discrete target space. Namely, all multiloop amplitudes (three or more loops) depend on 
the parameters of the boundaries through a product of factors associated with loops. In 
our case the differentiation in \i for fixed X±, X 2 gives 

d X] = J_J_ 1 r sinh[(g-l)^] S iDh[(p-l)g a ] 1 

^dfi XX{ U 21 dX 1 dX 2 cosYy{e 1 )-cosYv(6 2 ) V sinh(^i) smh(p6 2 ) J 

_ _d 8 sinh[(g - fc)6>i] sinh[(g - k)9 2 ] 

~ dXldX^f-^ sinh(^i) sinh(^ 2 ) 

(9.15) 

Similarly, the derivative of the YY correlator decomposes into a sum of p — 1 terms. 

Comparing ( |9.15|) with the standard expression for the two- loop correlator ( [37| , 



eq.(4.42)) we can express the operator P in the (p,q) critical regime of the two-matrix 
model as a linear combination of the order parameters of the corresponding string theory. 

25 



Indeed, the loop-loop correlator diagonalizes in the momentum space k/q,k = 1, 2, q—1, 
Introducing the loop operator ^k(X) with momentum kjq we have 

/jv/vw, /yu 9 9 sinh[(g - fe)gi] sinh[(g - k)6 2 } , . 

Mi al 2 sinh(gfc'i) smh(g02j 



This formula can be obtained from eq. (4.42) of ref.|37]], or, more easily, from the Laplace 
image of the expression for the three-loop correlator (eq. (4.58) of ref. ||37|l ) with one of 
the loops shrinked to a point. Comparing ( |9.15[ ) and ( |9.16| ) we can conclude that the 
Heisenberg momentum P can be written as a linear combination of the loop operators 
^^(X), k = 1, ...,q — 1. In the same way the operator Q is a linear of the loop operators 
\l/fc(y), k = 1, 2, — 1. Note that the XX ( YY) correlators depend only on q (p). This 
phenomenon occurs only for the two-loop correlator and is due to the Euler characteristic 
of the world-sheet being zero. The two-loop average, however, contains all the operators; 
this is related to the fact that the Euler characteristic of the cylinder is zero. As for the 
mixed correlator Wxy, its derivative with respect to A cannot be written in a factorized 
form in the non-symmetric case. 

All these formulae have an especially nice interpretation in the Z2-symmetric case of 
the two-matrix model U = V, describing the (m, m + 1) unitary models, which we consider 
in more detail in the next section. 



10. Symmetric (V = U) realization of the unitary (m, m + 1) models 

In this case we can use the equation of motion ( |2.17| ) . The formula for the multicritical 
potential is given by ( |4.4| ) if one takes there p = q = m. 

If one uses only the even potentials (polynomials in the even powers of the matrices) 
one has to take, instead of eq. ( [4.3|) , the following critical expression for X*(z): 

d (l-z 2 '\ m ~ 1 

r x.W = (io.i) 
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which is the only one satisfying all the conditions of m-criticality and the behaviour for 
z — > 0. The expression for the multicritical potential is the same as for the general 
symmetric case. 

Some examples of these potentials are given in the Appendix. Note that at least 
the even (m = 2k) unitary cases provide a non-perturbatively stable definition of the 
corresponding rational string theory in the double scaling limit, since the corresponding 
multicritical potential contains the positive highest power coefficient. As far as we know, 
this is the simplest possible stabilization of the double scaling limit of the unitary models 
of gravity. 

The symmetric case provides a direct statistical interpretation of the two matrix model 
in terms of the Ising spins on the dynamical random lattice. In this case the operator Y is 
the Wj 2 image of X and commutes with it. The loop- loop correlators found in the previous 
section correspond to the boundary conditions with fixed directions of the Ising spins: 
Wxx and Wyy have the same directions of the spin along the two boundaries, whereas 
Wxy has opposite directions. The /(/-derivative of the loop-loop correlators factorizes both 
for Wxx = Wyy and Wxy- 

,, d W (X Y\- 3 3 V" siiih(fcfli) sinh(fcfl 2 ) 
^t-W'm Ai, X 2 ) - > — -—r — —, — — (1U.2) 

oil oX\ oX 2 ^— J x smh(m6>i) smh(m^2j 

,* 3 W (X X\- 3 3 VV U sinh (^i) sinh(fcg 2 ) nn ,x 
o/i oXi oX 2 j—* sinh(mc'i) sinh(mt'2) 



where 



X 1 = ii m cosh(m6'i),X 2 = // m cosh(m0 2 ) (10.4) 



We see that the two-point correlator diagonalizes if we pass to the operators imposing 
symmetric (+) and antisymmetric (— ) boundary conditions. 
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11. Discussion 

We made an attempt to explicitly work out the idea proposed by M. Douglas about 
the possibility to represent all models of (p, q) rational matter interacting with 2D gravity 
as the multicritical regimes of the two matrix model, which describes random planar graphs 
of specially adjusted statistics with the Ising spins on it. Although it should be clear to 
the reader that this purpose has not been achieved completely in our paper, we made a 
few steps forward in this direction, leaving many other important details for the future. 

We determined the pairs of (p, q) multicritical potentials, i.e. the sets of coupling con- 
stants ( "temperatures" ) rendering the (p, q) matter critical. By shifting the cosmological 
coupling away from its critical value we reconstructed the known scaling solutions for the 
one loop averages in these models. 

We found an intriguing relation between the canonical momentum and the loop aver- 
age in the dispersionless (large N) limit of the generalized KdV formalism. 

The two-loop correlators are given in a closed form, which is valid for arbitrary poten- 
tials. Some of our formulae and the method of calculation apply for any chain of matrices 
admitting orthogonal polynomial representation. The scaling form of loop-loop correlators 
is compatible with the previous results obtained for the RSOS models coupled with gravity: 
in the unitary (m, m + 1) model the correlator can be factorized into (m — 1) terms describ- 
ing the corresponding finite number of excitations (corresponding to the order operators) 
propagating between the macroscopic states (loops). 

One of the most attractive features of the construction of the theory of {p, q) matter 
coupled with 2D gravity via the two matrix model, is the explicit (p, q) versus (q, p) duality. 
This last property is far from evident in the standard construction based on the RSOS 
models on a random lattice. 

Among the problems to be addressed further are the following: 

l.The calculation of the multi-loop correlators in form as explicit as was done for the 
one matrix model. 
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2. The description of the double scaling limit and of the topological 1/N expansion, 



using the exact (for any N) operator equations (|2.15|) . This becomes especially interesting 
in view of the possibility to give within this model a nonperturbative stable definition of 
the double scaling limit of (2n, 2n + 1) unitary models (note that the critical potentials 
in the Appendix (even case) have a positive coefficient in front of the highest power for 
m = 2n). 

3. It would be interesting to compare the 1/N expansion in the two matrix model 
with the "surface surgery" procedure worked out for the ADE models on random surfaces 



37f] . In particular, it would be very instructive to reproduce the open string amplitudes 
found recently in the framework of the SOS model P0| . m the two matrix model they 
must correspond to the loop amplitudes with the mixed boundary (including the products 
of both matrices) in the loop operator. Some recent results in this direction were obtained 



in 41] using the methods of ref. [42 



4. A possible application of the formalism worked out here could be the investigation 
of the critical regimes of the so-called induced gauge theory |2£] in the large N limit. 
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Appendix A. (p,q) Multicritical Potentials for the Two Matrix Model 

A.l. Nonsymmetric case 



q = 2,p = 3 



3,p 



q = 4, p = 5 



2,p = 4 



q = 3,p = 5 



q = 5,p = 7 



V(y 
U(x 
V(y 
U(x 
V(y 
U(x 
V(y 
U(x 
V(y 
U(x 
V(y 
U(x 
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= -7x + 
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— + T 



-23y-23y 2 + — ^ 
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4 5 

115 x 3 345 x 4 161 x 5 
-23 x H h 



23 x 6 x 7 



(A.l) 



3 4 ' 5 

A. 2. Symmetric case U = V = V m describing the (m, m + 1) -unitary models p = q = m 



V 3 (x) 



-3x — 



3x 2 x s 



2 + 3 



V 4 (x) = 8x + 2x 2 + 



x 

T 



V 5 (x) = -15 x 



5x< 



+ 15 x" 



15 x 4 



x'- 



V Q {x) = 24x + 3x 2 + 



4 

+ 39 x 4 + 



224 x 3 „„ , 24 x 5 



x 6 

5 + Y 



(A.2) 
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A. 3. Even Symmetric Potentials V m for (m, m + 1) models 



V s (x) = 2 x 2 - a 



12 

2 „4 



X 6 



V 4 (x) = ll^-x 4 + 

w , 452 4 4x 6 (A ' 3) 

V5{X) = 13 X + 



5 56 

, . 64475 x 2 604 x 4 158 5x 8 x 10 
^) = — 63 3" + ^ T + W 
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